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SUMMARY 


Recent  investigations  on  the  influence  of  edge  conditions  on  the  critical  axial  load  of 
cylindrical  shells  are  extended  here  to  include  cylinders  with  elastic  edge  restraint. 
It  is  shown  that  weak  edge  conditions  are  unlikely  to  affect  the  critical  load  signifi¬ 
cantly  in  experiments  or  in  practice.  Also  the  effects  of  symmetric  initial  imper¬ 
fections  were  investigated.  The  critical  load  for  cylinders  with  initial  imperfections 
in  the  form  of  a  cosine  function  was  found  to  be  somewhat  lower  than  was  indicated  in 
a  previous  analysis  by  Koiter.  The  analysis  also  disclosed  that  very  short  cylinders 
would  not  be  sensitive  to  the  types  of  imperfections  which  were  considered.  A  few 
tests  were  performed  and  their  results  tend  to  support  this  conclusion.  { 
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NOTATION 


F. 

1 

H 


L 

M 


integration  constants 

Et3/[  12(1  -  v2)} 

YoungTs  modulus 
stress  function  [See  Eq.  (21)] 
values  of  F  at  discrete  point 
radial  force  at  shell  edge  [See  Eq.  (6)] 
shell  length 

total  number  of  discrete  points 


Mx,  Mxy,  Myx  bending  and  twisting  moments  per  unit  width 
N  compressive  axial  load  per  unit  width 


N/(2yEt) 

normal  and  shear  forces  per  unit  width 

value  of  N  in  prebuclding  range 
x 

critical  axial  load  on  cylinder 
lateral  displacement  [See  Eq.  (21)] 


W.  value  of  W  at  discrete  point 

Z  L2(l  -  Al/2/(rt) 

ax  k(l  +  N)l/2 

—  1/2 

a,,  k(l  -  N)x/ 

b  length  of  initial  imperfection  wave  [See  Eq.  (34)] 
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b. 


d. 

1 


e. 

1 


f 

h 

i,  m,  q 

k 

n 

P 

P 

r 

t 

u,  v,  w 


w 


"1 
X,  y 


a 


a 


y 


coefficients  [See  Eq.  (IS)] 

2 

end  ring  area  divided  by  t 

coefficient  in  displacement  function  [See  Eq.  (29)] 
coefficients  [See  Eq.  (18)] 
coefficients  [See  Eq.  (20)] 
stress  function 

distance  between  neighboring  discrete  points 
integers 

<2yf1/2 

number  of  waves  in  circumferential  direction 

internal  pressure 

(pr)/(Et) 

shell  radius 

shell  thickness 

axial,  circumferential,  and  radial  displacements 
value  of  w  in  prebuckling  range 
initial  displacement 

value  of  initial  displacement  [See  Eq.  (32)] 
axial  and  circumferential  coordinates 
value  of  F  or  W  [See  Eq.  (23)] 

(aiL)/(2r) 

(a2L)/(2r) 

coefficient  in  prebuckling'  displacement  [See  Eq.  (14)] 
(t/r)/[  12(1  -  v2)\l/2 

viii 


A 

length  of  initial  imperfection  wave  [See  Eq.  (32)] 

iir  r/L 

\ 

spring  constants  [See  Eq.  (7)] 

Xl’  X3 

nondimens ional  spring  constants  [See  Eq.  (8)] 

X. 

l/X 

1 

1 

Mi 

nondimens ional  amplitude  of  initial  displacement  component 

M 

nondimens  ional  amplitude  of  sinusoidal  initial  displacement 

V 

Poisson’s  ratio 

Cl) 

initial  displacement  [See  Eqs.  (33,  34)] 

aCR 

critical  axial  stress 

ctCL 

classical  value  of 

V4 

34/ 3x4  +  234/3x2  3y2  +  34/3y4 

When  subscripts 

x  and  y  follow  a  comma,  they  indicate  partial  differentiation  of  the 

principal  variable  with  respect  to  x  or  y  .  Primes  indicate  total  derivatives  with 

respect  to  x  . 

ix 

INTRODUCTION 


Due  to  the  lack  of  agreement  between  theory  and  test,  empirical  curves,  based  on  a 
large  number  of  tests,  are  usually  employed  in  the  design  of  axially  compressed  cyl¬ 
inders.  However,  many  of  the  test  specimens  were  probably  manufactured  poorly 
and  the  design  curves  are  generally  conservative.  Of  course,  as  more  test  results 
are  made  available,  questionable  test  points  may  be  disregarded  and  better  design 
curves  may  be  constructed.  However,  the  designer  will  still  be  on  uncertain  ground 
until  the  buckling  problem  for  the  axially  loaded  cylinder  is  fully  understood.  The 
effects  of  initial  imperfections  and  of  differences  in  edge  conditions,  for  instance, 
cannot  be  satisfactorily  evaluated  through  experimental  analysis  alone. 

The  influence  of  edge  conditions  as  a  possible  explanation  of  the  discrepancy  between 
theory  and  test  was  explored  by  Nachbar  and  Hoff  (Ref.  1),  Stein  (Ref.  2),  and 
Fischer  (Ref.  3).  Nachbar  and  Hoff  considered  the  case  of  free  edges  and  found  the 
critical  load  to  be  37  percent  of  the  classical  buckling  load.  Stein  presented  an  analysis 
of  cylinders  with  supported  edges  which,  for  the  first  time,  included  an  accurate  pre¬ 
buckling  analysis.  For  the  particular  set  of  boundary  conditions  considered  by  Stein, 
the  computed  critical  load  was  about  half  the  classical  buckling  load.  Fischer  pre¬ 
sented  a  similar  analysis  but  with  different  in-plane  boundary  conditions.  Fischer's 
critical  load  was  as  high  as  85  percent  of  the  classical  buckling  load.  With  boundary 
conditions  identical  to  those  by  Stein  but  with  a  membrane  prebuckling  solution  Ohira 
(Ref.  4),  Hoff  (Ref.  5),  and  Hoff  and  Rehfield  (Ref.  6)  obtained  results  close  to  those 
by  Stein.  The  investigations  by  Stein  and  by  Fischer  were  extended  by  Almroth  (Ref.  7) 
to  cover  eight  different  sets  of  boundary  conditions.  It  was  shown  that,  except  in  the 
free  edge  case,  a  substantial  reduction  of  the  buckling  load  occurs  only  when  the  con¬ 
ditions  of  simple  support  and  of  zero  tangential  displacement  are  applied.  Neither 
the  case  of  zero  transverse  shear  force  (Nachbar  and  Hoff)  or  zero  tangential  shear 
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force  (Stein)  is  likely  to  be  completely  realized  in  tests  or  in  practical  applications, 
but  it  is  possible  that  considerable  reduction  may  be  due  to  elastic  restraint.  There¬ 
fore,  a  buckling  analysis  with  elastically  supported  edges  is  presented  here. 

Although  weak  edge  restraint  and  also  disturbances  such  as  accidental  lateral  loads 
may  be  of  importance  in  some  practical  applications,  these  factors  are  certainly  of 
minor  importance  in  tests  performed  under  laboratory  conditions.  Consequently, 
the  presence  of  initial  geometrical  imperfections  still  appears  to  be  the  major  reason 
for  the  discrepancy  between  theory  and  test. 

Effects  of  imperfections  in  the  cylinder  were  investigated  by  Donnell  and  Wan  (Ref.  8). 
Their  nonlinear  analysis  demonstrates  clearly  the  importance  of  the  initial  imperfec¬ 
tions,  but,  due  to  the  approximations  involved,  it  does  not  appear  possible  to  use 
their  analysis  directly  for  design  purposes. 

Some  analyses  consider  cylinders  with  axially  symmetrical  initial  imperfections  only. 
Koiter  (Ref.  9)  first  presented  a  general  theory  which  was  valid  only  for  imperfections 
of  small  amplitude,  and  later  a  less  restricted  theory  (Ref.  10).  Babcock  and  Sechler 
(Ref.  11)  derived  a  solution  to  the  problem  which  is  quite  similar  to  Koiter's  latest 
theory.  Both  analyses  consider  only  the  case  in  which  the  initial  displacements  are 
described  as  a  simple  trigonometric  function  of  the  axial  coordinate.  Babcock  and 
Sechler  restricted  their  analysis  to  the  case  in  which  the  axial  halfwave  length  of  the 
initial  displacement  equals  the  shell  length.  Koiter  restricted  his  numerical  analysis 
to  include  only  imperfections  which  are  proportional  to  the  symmetrical  buckling 
pattern  of  an  infinitely  long  perfect  cylinder.  Also  Fischer  (Ref.  3)  considered  the 
effect  of  axially  symmetrical  initial  imperfections .  The  approach  used  by  Fischer 
allows  more  freedom  in  the  choice  of  the  shape  of  the  imperfection.  However, 
numerical  results  are  presented  only  for  one  case. 

Unless  the  analysis  is  restricted  to  symmetric  initial  imperfections,  the  numerical 
work  will  require  the  solution  of  nonlinear  differential  equations.  Therefore,  it  is 
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reasonable  to  explore  first  the  symmetrical  case  fully  although  such  an  analysis  may 
not  constitute  a  final  solution  of  the  problem.  In  the  present  analysis,  only  axially 
symmetrical  imperfections  are  considered.  The  imperfections  are  represented  by 
trigonometric  series,  and  the  effect  on  the  buckling  load  of  imperfection  of  various 
shapes  is  studied. 
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BASIC  EQUATIONS 


The  so-called  Donnell  equations  are  used  in  the  present  analysis.  Thus  for  a  cylinder 
with  an  initial  imperfection,  w  , 


N  = 


x 


N  = 


\(  1  2  1 

/  w 

12—  \] 

(u  +  -  W  +  w  w  J 

1  +  Wv  -  —  + 

z  w  +  w  w  II 

IV  ,x  2  ,x  ,x  ,x/ 

V  ,y  r 

2  >y  .y  .y/J 

'  7^7  Ku,x  +  Iw2x  +  w,xw,x)+  (v, 


w  1  2 

+xw  +  w  .. 
y  r  2  ,y  , y  ,y> 


w  )| 

»y  »y/i 


N 

xy 

M  = 
x 

M  = 


Et 


M 


xy 


2(1  +  ,)  I“,y  +  v,x  +  \yw,x  +  w,xw,y  + 


-  D(w  +  vw  ) 
,xx  ,yy' 


-  D(^w  +  w  ) 
v  , xx  , yy7 


M  =  D(1  -  v)  w 
yx  ,xy 


(1) 


By  use  of  these  equations  the  total  potential  energy  of  the  loaded  cylinder  can  be 
expressed  in  terms  of  the  displacement  components  u  ,  v  ,  and  w  .  From  this  energy 
expression  we  find  through  the  variational  approach  the  three  equilibrium  equations: 


N  +  N 
x,x  xy,y 


N  +  N 
y,y  xy,x 


=  0 


=  0 


DV4W 


N(w  +  w  )  +  N  (w  +  w  ) 
x'  j  xx  ,xx'  yv  ,yy  ,yy; 


+  2N  (w  +  w  ) - N+p 

xy  ,xy  ,xy7  r  y  * 


(2) 
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An  Airy's  stress  function  f  is  introduced  such  that 


Nx  ■  f,yy  ; 


N  =  f  ; 
y  .xx 


N  -  -  f 
xy 


xy 


(3) 


The  first  two  equations  may  then  be  replaced  by  the  compatibility  equation 


—  v^f  =  -  w  w  +  2w  w  -w  w  -w  w  +  —  w  (4) 

EtV1  W, xy  w,xxw,yy  ,xyw,xy  ,xx,yy  ,yy  ,xx  r  ,xx 

Since  the  analysis  is  restricted  to  axially  symmetrical  initial  imperfections,  the 
equilibrium  and  compatibility  equations  may  be  written  as 


DV4w  =  f5yy(wjXX  +  WjXX)  +  f jXxw,yy  -  2f,xyw,xy  ■  Ff.xx  +  P 

—  v4f  =  w2  -w  w  -  w  w  +  i  w 
Et  ,xy  ,xx  ,yy  ,xx  ,yy  r  ,xx 


(5) 


In  addition  to  the  equilibrium  equations,  the  variational  approach  yields  the  following 
natural  boundary  conditions  at  x  =  ±  L/2: 


=  N  =  -  N 

.yy  x 

=  n  =  o 

.xy  xy 

Mx  =  -D<w,xx+  ^.yy)  =  0 
H  *  f,yy<w,x  +  »,x>  -  f,xyw,y  '  r"  lw,xxx  +  <2  '  ”>  w,xyy> 


>■  (6) 


=  0 
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All  admissible  sets  of  boundary  conditions  are  given  by 


w 


w 

,x 

v 

u 


=  X11H  +  X12Mx  +  X13Nxy  +  X14Nx 

X21H  +  X22Mx  +  X23Nxy  +  X24Nx 
X31H  +  ^32Mx  +  A33Nxy  +  X34Nx 

=  X41H  +  X42Mx  +  X43Nxy  +  X44Nx 


(7) 


It  will  be  sufficient  for  the  purpose  of  this  analysis  to  consider  only  the  effects  of  the 
diagonal  terms  in  Eqs.  (7).  Consequently  in  the  numerical  analysis  we  assume 
A  =  0  if  i  ^  j  .  Also,  as  the  choice  of  boundary  condition  for  loads  or  displace¬ 
ments  in  the  axial  direction  generally  has  little  effect  on  the  critical  load  we  will  use 

N  =  -  N  at  x  =  ±  L/2 
x 


Nondimentional  spring  constants  are  defined  as  follows: 


A^  =  Et/r 


X3  X33  Et^r 


Thus  the  boundary  conditions  may  be  written  (for  total  displacements) 


N  =  -  N 
x 

w  +  A  (Hr)/(Et)  =  0 

wx  +  A2Mx/(Et2)  =  0 
v  -  A3(Nxyr)/Et  =  0 


(8) 


(9) 
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PREBUCKLING  DISPLACEMENTS 


As  only  axially  symmetric  imperfections  are  included  in  the  analysis  the  prebuckling 
displacements  will  be  independent  of  the  circumferential  coordinate.  Derivatives  of 
prebuckling  displacements  with  respect  to  y  will  vanish  and  from  Eqs.  (5)  we  obtain 


^O.kxxx^Vxx'^'Vxx-P  =  -Nw 


(l/Et)fQ  xxxx  (i/i)  wojXX 


where  the  subscript  zero  is  used  to  signify  prebuckling  quantities,  w  represents  the 
initial  lateral  displacement,  and  N  -  -  Nx  . 

Prebuckling  displacements  must  be  symmetric  about  x  =  0  ,  and  therefore  only  even 
functions  are  included  in  the  solution  of  the  differential  equations  and  the  boundary 
conditions  need  only  be  enforced  at  x  =  +  L/2  .  The  homogeneous  solution  consequently 

is 
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and  A j  ,  are  constants  of  integration.  It  is  assumed  that  the  initial  displace¬ 
ments  may  be  expressed  as  a  trigonometric  series  as  follows 


(13) 


Consequently,  Eqs. 


(10)  have  a  particular  solution 


where 


2NyM  (t/r)  f? 

~  —  2  2  4 

1  -  2Nyff  +  y 

, .  i7rr 

^i  L 

y 


(14) 


Due  to  the  existing  symmetry  only  two  boundary  conditions  are  needed,  the  second  and 
third  of  Eqs.  (9).  For  the  case  of  axial  symmetry  they  may  be  written 


wA  -  rX,  [2yN(wA  +  w  )  +  r2y2wA  1  =  0 
0  1  l  v  0,x  ,x'  '  0,xxxj 


wA  -  rX0wn  ~  0 
0,x  2  0,xx 


(15) 


By  addition  of  the  previously  found  solution,  Eq.  (11),  and  the  particular  solution, 
Eq.  (14),  and  substitution  of  the  sum  into  the  boundary  conditions,  Eq.  (15),  we  find 
for  the  prebuckling  displacements 


w 


0 


(- 


Et  Et/ 


1  +  Aj  sin 


+  A 


2 


(16) 
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where 


e2d2  d4  ^d5  +  el^ 
d2d3  '  dld4 


d2d3  -  dld4 


bl  -  Xl(*\  +  2yNb3) 

b2  -  Xl(At  + 


b3  "  \J  b! 


b4  "  X27  (r//t)  b( 


=  -  (2vyN  +  p) 


=  (aiL)/(2r) 


(a2L)/(2r) 


sin  a  ^  sinh  a2 


cos  a1  cosh  a  ^ 


|(1  -  N)^*2  sin  cosh  a2  +  (1  +  N)  ^  cos  sinh  o' 2 j 
k  |(1  +  N)1^2  sin  cosh  o;2  -  (1  -  N)1^2  cos  a  1  sinh  a  2 
2k2  [n  sin  sinh  «2  -  (1  -  N2)1//2  cos  a1  cosh  a 2] 


—  2  1/2  1 
sinh  a  2  -  (1  -  N  )  cosaj^  cosh  a2\ 


-  2k2  |(1  -  n2)1^2  sin  a  ^  sinh  a  1  +  N  cos  a  x  cosh  a 2 


[(1  -  N)1/2  (1  +  2N)  sin  a  ^  cosh  a2  -  (1  +  N)1/2  (1  -  2N)  cos  ax  sinh  aj 


3  (1  +  N)1,/2  (1  -  2N)  sin  a1  cosh  a2  -  (1  -  N)1//2  (1  +  2N)  cos  oij  sinh  a2J 


9 


[  wo/r  +  xi  [ 


r  y  w„  +  2yNw»  -  2yNw 

0,xxx  0,x  ,x 


_  r  S  [^4i+2  ^4i  +  Xly  (f4i+l/34i+l  f4i 


4i+3K4i+3, 


+  X,  2yN 


{?4i+l  [^41+1  (t/r)  ^41+1]  ^4i+3  ^4i+3  (t/r)  ^41+3!}] 


[-W0,x  +  rVr/t)  ^O.xx. 


x  =  L/2 


r  ^  [?4i+3^4i+3  ”  ^4i+1^4i+l  +  A2(r/^  y  (^4i^4i  ^4i+2^4i+2)J 
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BUCKLING  ANALYSIS 


Buckling  equations  are  derived  from  the  nonlinear  differential  equations  through 
decomposition  of  stresses  and  displacements  into  prebuckling  and  incremental 
quantities.  Then  in  the  resulting  equations,  the  terms  containing  only  prebuckling 
quantities  may  be  subtracted  out  by  virtue  of  prebuckling  equilibrium,  and  non¬ 
linear  terms  in  the  (infinitesimal)  incremental  quantities  may  be  disregarded. 

This  procedure  results  in  the  replacement  of  nonlinear  equations  of  equilibrium 
by  linear  equations  of  stability. 

The  space  variables  in  these  equations  may  be  separated  by  means  of  the  following 
representation  for  f  and  w  : 


f  =  F  cos  (ny/r) 
w  =  W  cos  (ny/r) 


(21) 


The  problem  is  thus  reduced  to  the  solution  of  a  set  of  linear  ordinary  differential 
equations  with  variable  coefficients.  A  finite  difference  approach  will  be  used  to 
effect  a  numerical  solution.  The  differential  equations  and  the  boundary  conditions 
are  converted  into  a  system  of  algebraic  equations  by  use  of  the  following 
appr oxim  ations : 


y’  =  [yi+1  - 

y"  =  Iyi+1  -  2y,+  y^l/h2 
y!"  =  [i/2  y1+2  -  y1+1  +  y,.!  -  1/2  yi-2I/h3 
=  iyi+2  -  4Vi 4  %  -  “/i-i +  V21/h4 
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The  unknowns  of  the  problem  are  the  values  of  \Y  and  F  (W .  and  Fp  at  a  finite 
number  of  points,  including  the  fictitious  points  outside  the  cylinder  which  are  needed 
for  definition  of  the  boundary  conditions.  The  points  are  spaced  at  a  constant  dis¬ 
tance  h  .  The  total  number  of  points  is  M  =  L/(2h)  +  5  . 

The  number  of  unknowns  equals  2M  .  The  equilibrium  and  compatibility  equations 
at  each  of  the  real  points  on  the  cylinder  and  the  eight  boundary  conditions  provide  2M 
equations.  For  convenience  in  the  numerical  analysis,  the  unknowns  (W\  and  F '.) 
are  denoted  by  y.  such  that 


y0  =  w 

J2m  m 

(23) 

^2m-l  ^m 

The  finite  difference  equations  of  equilibrium  and  compatibility  are  identical  to 
Eqs.  (31)  and  (32)  of  Ref.  7,  except  that  (WJJ)  is  replaced  by  (WJJ)  +  . 

Due  to  symmetry  of  loading  and  boundary  conditions,  the  incremental  displacements 
are  either  symmetric  or  antisymmetric  about  x  =  0  .  For  the  symmetrical  case, 
the  boundary  conditions  at  mid-length  are  given  by  Eqs.  (30)  in  Ref.  7.  In  the  case 
of  antisymmetry  the  boundary  conditions  at  x  =  0  are 

W  =  0 
W"  =  0 

(24) 

F"  =  0 
v  =  0 

It  can  be  shown  that  the  condition  v  =  0  is  identical  to  the  condition  that 
2 

F"  +  vn  F  -  0  .  By  virtue  of  the  third  condition  of  Eqs.  (24),  then  v  =  0  can 
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here  by  replaced  by  F  =  0  .  In  terms  of  finite  differences  the  boundary  conditions 
for  the  antisymmetric  case  are 


W3  =  *6  =  ° 

(W")3  =  (y8  -  2y6  +  y4)/h2  =  o 
(F")3  =  (y7  -  2y5  +  y3)A2  =  0 


(25) 


The  finite  difference  expressions  for  the  elastic  boundary  conditions  at  the  shell 


edge  (Eqs.  9)  in  terms  of  finite  differences  are 
<l/n2>  Nx  =  y2M-5  = 0 

W  +  lfeXl  H  =  y2M-4  “  Al|r  [2y^  +  y2(2  "  ^)](y2M-2  "  y2M-6)/(2h) 

+  y2r2(l  y2M  “  y2M-2  +  y2M-4  “  2  y2M-e)/h  }  ~  ° 

W,X  +  Et2  *2  ^  =  (Y2M_2  ”  y2M-6)/,(2h) 

o 

,  2\  „  .  /,  2  vn 

-  X2(r/t)  y  |r(y2M_2  “  2^2M-4  +  y2M-6Vh  r  y2M-4 

v  “  >3  Nxy  =  y2M-4  "  r  <y2M-3  '  2y2M-5  +  y2M-7^h 

-  ^  y2M-5  +  A3  rn2(y2M-3  ‘  y2M-7)/2h  =  ° 


The  equilibrium  and  compatibility  equations  together  with  these  boundary  conditions 
define  a  linear  homogeneous  equation  system  of  size  2M.  This  equation  system  can 
have  nontrivial  solutions  only  if  the  determinant  of  the  coefficients  equals  zero.  The 
lowest  value  of  N  for  which  the  determinant  vanishes  corresponds  to  the  critical  load 
of  the  cylinder.  This  value  of  N  will  be  found  by  a  trial  and  error  method. 
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RESULTS 


Perfect  cylinders 

Numerical  results  were  obtained  from  an  IBM  7094  digital  computer.  The  accuracy 
of  these  results  depends  on  the  use  of  a  sufficiently  fine  mesh  size  in  the  finite  differ¬ 
ence  scheme.  The  investigations  of  Ref.  7  have  essentially  established  the  number  of 
points  needed  for  satisfactory  accuracy  and  therefore  the  number  of  points  was  varied 
here  only  in  a  few  check  cases. 

The  limiting  cases  of  complete  restraint  are  obtained  if  the  corresponding  non-dimen¬ 
sional  spring  constants  are  set  equal  to  zero.  In  order  that  free  edge  conditions  could 
also  be  included  in  the  analysis,  the  program  was  branched  such  that  the  inverse  of 
these  constants  could  be  used.  Complete  agreement  was  found  with  the  results  of 
Ref.  7  for  perfect  cylinders  in  the  limiting  cases. 

Figure  1*  shows  the  variation  of  the  critical  load  with  shell  length  for  the  case  in  which 
the  edges  are  free  in  the  radial  direction.  This  case  was  not  included  in  Ref.  7. 
Results  for  long  shells  are  available  in  Ref.  1  for  this  edge  condition.  Although  the 
boundary  conditions  in  the  two  analyses  were  somewhat  different,  the  numerical  re¬ 
sults  are  in  good  agreement. 

Buckling  loads  for  shells  with  elastic  edge  restraint  are  shown  in  Figs.  2,  3,  and  4. 
The  variation  of  the  critical  load  with  the  tangential  restraint  parameter  A  is  shown 

O 

in  Fig.  2  for  different  values  of  the  number  of  circumferential  waves.  Curves  which 
fall  above  their  neighbors  in  the  entire  range  are  generally  not  shown  in  order  that 
the  figure  be  as  clear  as  possible.  When  the  rotational  edge  restraint  A 2  was  varied, 
it  was  assumed  that  the  shell  was  free  in  the  tangential  direction  as  otherwise  the 
spread  of  the  values  of  the  critical  load  over  the  Arrange  would  be  quite  small.  The 
results  are  shown  in  Fig.  3.  The  variation  of  the  critical  load  with  radial  edge 


^Figures  are  presented  at  the  end  of  this  section. 


restraint  X1  is  shown  in  Fig.  4.  Only  cylinders  with  v  =  0  are  considered  because 
with  N  =  0  at  the  edges,  the  critical  load  is  almost  independent  of  the  amount  of 
elastic  restraint.  From  these  figures  it  is  seen  that  the  curves  for  different  numbers 
of  waves  are  sometimes  arranged  in  a  disorderly  manner.  The  reason  for  this  is  that 
the  critical  load  as  a  function  of  the  wave  number  may  have  two  minima  as  illustrated 
in  Fig.  5. 

The  results  shown  in  Figs.  2,  3,  and  4  can  be  used  to  determine  the  critical  load  of  a 
perfect  cylinder  whenever  the  stiffness  of  the  substructure  is  known.  In  order  to 
facilitate  a  quantitative  understanding  of  the  curves,  buckling  loads  have  been  computed 
for  cylinders  with  an  end  ring.  Two  different  cases  have  been  considered.  In  the 
first  case,  the  cylinder  edge  is  supported  in  the  radial  direction  and  the  torsional 
stiffness  of  the  end  ring  is  omitted  such  that  it  provides  support  only  against  tangential 
displacement.  In  the  second  case,  the  torsional  stiffness  of  the  ring  is  included  but 
the  radial  support  is  removed. 

The  end  ring  cross  section  is  assumed  to  be  square.  For  such  a  ring  with  area  equal 
o 

to  ct  ,  the  nondimensional  elastic  constants  are 


*1 

=  l/Xj 

X2 

-  1A2 

0.  141  2  2  /t\2 

“  2(1  +  v)  C  n  \r)  j 

► 

(27) 

X3 

=  iA3 

-  "2 ■(*) 

There  is,  of  course,  also  a  coupling  term  between  tangential  and  radial  displacements 
but,  according  to  previous  assumptions,  this  term  is  omitted. 


The  variation  of  the  critical  load  with  the  size  of  the  end  ring  is  shown  in  Figs.  6  and 
7.  For  the  case  in  which  the  ring  is  radially  supported  (Fig.  6)  as  well  as  the  case 
in  which  the  ring  is  free  (Fig.  7),  it  is  seen  that  a  very  small  ring  is  sufficient  to 
make  the  critical  load  of  the  shell  appproximately  equal  to  the  critical  load  for  a  shell 
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with  complete  restraint.  In  both  cases  the  restraint  against  axial  displacements  was 
omitted. 

Figures  8,  9,  and  10  show  buckling  modes  for  cylinders  with  different  boundary 
conditions. 

The  mode  shape  shown  in  Fig.  8  corresponds  to  a  simply  supported  cylinder  v  =  0. 
Parameter  values  are  the  same  as  those  considered  by  Fischer  in  Ref.  3.  Although 
the  critical  loads  found  in  the  two  analyses  are  identical,  there  is  a  slight  difference 
in  the  mode  shapes,  inasmuch  as  the  corresponding  curve  in  Ref.  3  appears  to  have  a 
horizontal  tangent  at  the  end  of  the  cylinder.  Figure  9  shows  the  mode  shape  for  a 
cylinder  with  the  edges  free  in  the  tangential  direction  only.  In  Fig.  10  is  demon¬ 
strated  how  the  mode  shapes  vary  with  the  size  of  an  end  ring.  The  curve  with 
c  =  0  shows  the  mode  shape  for  a  cylinder  with  free  edges  and  the  curve  for  c  =  16 
approximates  the  mode  shape  for  a  clamped  cylinder. 

Some  consideration  was  given  to  the  question  of  which  boundary  conditions  are  generally 
applicable  in  practice.  A  large  number  of  cylindrical  shells  constituting  structural 
parts  of  missiles  and  spacecraft  were  studied.  In  view  of  the  above  results  for  cylin¬ 
ders  with  end  rings,  it  was  concluded  that  in  no  case  would  the  critical  load  of  a 
perfect  cylinder  be  significantly  below  the  critical  load  for  the  simply  supported 
cylinder  with  v  =  0  and  w  =  0.  It  is  judged  that  cylinders  with  extremely  weak  edge 
conditions  will  rarely  occur  in  practice. 

Cylinders  with  imperfections 

Koiter  in  Ref.  9  gives  the  critical  load  of  infinitely  long  cylinders  with  an  initial  dis¬ 
placement  in  the  form  of  a  simple  trigonometric  function.  According  to  this  analysis, 
the  critical  number  of  waves  is  given  by 

n2  =  \  Vs(l-  r2)  (r/t)  ,  (p< 0.3)  (28) 
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For  reasonably  long,  simply  supported  shells  (N  =  0  ,  v  =  0)  ,  and  with  this  value 
of  n  ,  it  was  found  that  the  numerical  results  from  the  present  analysis  are  in  good 
agreement  with  those  by  Koiter.  However,  in  the  present  analysis,  lower  critical 
loads  were  found  for  other  values  of  n  .  For  instance,  in  the  case  of  a  cylinder  with 
r/t  -  100  and  the  amplitude  of  the  imperfection  given  by  /x  =  0.  2,  Koiter's  analysis 
as  well  as  the  present  analysis  indicates  a  critical  load  equal  to  0.447  times  the 
classical  buckling  load  when  n  =  9  .  The  lowest  buckling  load  which  occurred  at 
n  =  4  was  found  to  be  .389  times  the  classical  buckling  load.  The  mode  shapes 
corresponding  to  n  =  4  and  n  =  9  are  shown  in  Fig.  11.  It  is  interesting  to  note 
that  with  n  =  9  ,  the  cylinder  has  nodes  at  all  points  of  maximum  outward  initial  dis¬ 
placement.  Koiter  assumed  that  the  buckled  shape  could  be  represented  by 

oo 

w(x,y)  =  t  c.  cos  |(2i  -  1)  cos  (29) 

i  =  1 

for  the  case  in  which  the  initial  displacement  is  given  by 

w(x)  =  -  nt  cos  (^~J  (30) 

This  buckling  mode  is  in  good  agreement  with  the  mode  found  here  for  n  =  9  but 
evidently  Koiter's  analysis  is  too  restrictive  to  include  the  buckling  pattern  corres¬ 
ponding  to  n  =  4  .  Therefore,  Koiter's  analysis  was  modified  here,  such  that  the 
incremental  displacements  are  now  represented  by 

OO 

w(x,y)  =  t  ^  c.  cos  [(2i  -  1)  ^  cos  (31) 

i  =  1 

Koiter1  s  original  analysis  (m  =  1)  assumes  that  nodes  in  the  lateral  displacement 

occur  at  each  point  of  maximum  outward  initial  displacement.  With  the  modified 

th 

displacement  function,  the  nodes  occur  at  every  m  maximum  point.  The  modified 
formulation  was  programmed  for  the  digital  computer  and  m  was  varied  to  minimize 
the  critical  load.  The  results  so  obtained  were  in  good  agreement  with  results  from 
the  finite  difference  analysis  for  reasonably  long  shells .  The  program  was  used  also 


to  determine  the  critical  load  for  long  cylinders  as  a  function  of  the  amplitude  of  the 
initial  imperfections.  The  results,  which  are  valid  for  all  r/t  values,  are  shown  by 
the  middle  curve  in  Fig.  12.  It  was  found  for  p  =  0.2  that  the  critical  load  has  a 
minimum  at  m  =9.  This  means,  of  course,  that  the  critical  axial  wavelength  is 
finite  and  that,  for  sufficiently  long  shells,  the  finite  difference  approach,  also  with 
n  =  4,  should  lead  to  a  nodal  buckling  pattern.  This  was  shown  to  be  the  case  for  a 
shell  with  L/r  =  8.4,  and  the  corresponding  buckling  mode  is  shown  in  Fig.  13. 

Figure  14  shows  some  results  which  were  obtained  by  use  of  a  two-term  representa¬ 
tion  of  the  initial  displacements.  It  can  be  seen  that  with  a  fixed  amplitude,  the  buckl¬ 
ing  load  is  smaller  for  a  more  blunt  shape  of  the  initial  imperfection.  It  thus  appears 
reasonable  to  assume  that  a  lower  bound  to  the  critical  load  for  cylinders  with  periodic 
axially  symmetric  imperfections  of  fixed  amplitude  w  can  be  obtained  through  analysis 
of  the  case  in  which 

iA<x-^<(i+l)A  i  =  odd 

(32) 

iA  <  x  -  ^  <  (i  +  l)A  i  =  even 

Obviously  this  discontinuous  function  does  not  represent  meaningful  initial  displace¬ 
ments.  However,  a  finite  segment  of  the  Fourier  series  representation  of  this  function 
is  continuous.  With  increasing  number  of  terms  in  such  a  series  the  corresponding 
critical  load  increases  monotonically  but  within  a  narrow  load  range.  Consequently 
the  critical  load  based  on  a  one-term  series  may  be  considered  a  reasonably  close 
lower  bound  to  the  buckling  load  of  cylinders  with  periodic  initial  imperfections  of 
amplitude  wr  The  amplitude  of  the  first  term  is  about  1.25  w  and  the  lower  bound 
shown  in  Fig.  12  can  easily  be  constructed. 

For  long  cylinders  with  sinusoidally  varying  initial  imperfections  results  were  obtained 
also  with  other  boundary  conditions.  Figure  15  shows  the  critical  stress  versus  ampli¬ 
tude  of  imperfections  for  cylinders  with  the  edges  free  in  the  tangential  direction. 

Figure  16  shows  the  corresponding  stresses  for  cylinders  with  completely  free  edges. 
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It  may  be  seen  that  the  influence  of  edge  conditions  on  the  critical  stress  diminishes 
with  growing  amplitude  of  initial  imperfections . 

In  Ref.  3  Fischer  gave  a  solution  for  an  initial  imperfection  of  the  form 


w 


O)  COS^  (7rx)/(2b) 
0 


0  <  x  <  b 
b  <  x  s  L/2 


(33) 


Fischer's  numerical  results  indicate  that  an  initial  imperfection  as  described  by 
Eq.  (33)  is  considerably  less  critical  than  a  periodic  imperfection  of  the  same  amplitude. 
However,  the  numerical  analysis  was  very  restricted.  Only  two  values  of  b  were 
considered  and  only  positive  values  of  (o  were  considered,  and  it  appears  likely  that 
negative  values  (inward  bulge)  should  be  more  critical. 


Therefore,  the  influence  of  nonperiodic  imperfections  is  investigated  here.  The  form 
of  the  imperfection  was  chosen 


LO  COS^  (7rx)/(2b)  0  —  X  2= 

w  =  o  b  ^  x  < 

It  was  found  that  an  imperfection  of  this  type, for  all  values  of  b  ,  is  slightly  less 
harmful  than  a  periodic  imperfection  with  the  amplitude  U)/2  . 

It  appears  feasible  that  the  lower  curve  in  Fig.  12  really  represents  a  lower  bound 
even  if  cylinders  with,  nonsymmetric  imperfections  are  included.  If  this  could  be 
shown  to  be  the  case  by  use  of  suitable  experiments  or  by  additional  theoretical  analy¬ 
sis,  then  the  present  results  could  be  made  directly  useful  for  the  designer. 


b 

(34) 

L/2 
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Figure  17  shows  the  influence  of  initial  imperfections  on  the  critical  load  of  cylinders 
of  different  lengths.  It  appears  that  a  sufficiently  short  cylinder  is  not  sensitive  to 
initial  imperfections.  A  few  tests  of  ring  stiffened  cylinders  have  been  performed 
and  tend  to  confirm  this  conclusion. 

The  cylinders  were  manufactured  from  0.02  inch  aluminum  sheet  (AL  6061-T8)  and 
rectangular  aluminum  rings  were  bonded  to  the  skin  as  shown  in  Fig.  18.  Three 
cylinders  had  outside  rings  spaced  at  1.33  inches  and  three  cylinders  inside  rings 
spaced  at  1.  14  inches.  Three  cylinders  without  stiffening  were  tested  for  comparison. 
As  the  ring-stiffened  cylinders  are  such  that  buckling  between  rings  is  critical  the 
classical  buckling  load  for  all  three  different  types  of  cylinders  should  be 

PCR  =  2w  [3(1-  Arl/2  Et2  (35) 

Test  results  are  shown  in  Table  I.  The  values  of  Z  given  in  the  table  are  based  on 
the  ring  spacing  rather  than  on  the  shell  length.  It  can  be  seen  from  the  test  results 
that  the  ring- stiffened  cylinders  carry  the  classical  buckling  load  for  monocoque 
cylinders,  hence  also  the  experiment  indicates  that  very  short  cylinders  are  not  sen¬ 
sitive  to  initial  imperfections. 
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Critical  Loads  of  Cylinders  With  Elastic  Tangential  Restraint 
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ig.  2c  Critical  Loads  of  Cylinders  With  Elastic  Tangential  Restraint:  r/t  =  1000, 
L/r  =  0.222 
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Fig.  2d  Critical  Loads  of  Cylinders  With  Elastic  Tangential  Restraint:  r/t  =  1000, 
L/r  =0,76 
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Fig.  3a  Critical  Loads  for  Cylinders  With  Elastic  Rotational  Edge  Restraint: 
r/t  =  100,  L/r  =  0.7 


26 


CR  CL 


Fig.  3c  Critical  Loads  for  Cylinders  With  Elastic  Rotational  Edge  Restraint: 
r/t  =  1000,  L/r  =  0.222 
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Fig.  8  Buckling  Mode  for  Simply  Supported  Cylinder 


Fig.  11  Buckling  Modes  for  Imperfect  Cylinders 


Fig.  12  Buckling  of  Long  Cylinders  With  Periodic  Initial  Imperfections 
of  Amplitude  jut 


Fig.  13  Buckling  Mode  for  Long  Cylinder 


14  Influence  of  Shape  of  Initial  Imperfection  on  Critical  Load 
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Table  I 

TEST  OF  RING-STIFFENED  CYLINDERS 


Cyl.  No. 

Type 

B 

t 

Z 

p 

CR 

Theory 

Test 

1 

Monocoque 

15.96 

0. 0200 

1500 

16000 

10500 

2 

Monocoque 

15.76 

.0200 

1500 

16000 

12300 

3 

Monocoque 

15.71 

.0200 

1500 

16000 

8900 

4 

Inside  rings 

15.79 

.0200 

5.2 

16000 

15300 

5 

Inside  rings 

15.76 

.0205 

5.2 

16800 

14880 

6 

Inside  rings 

15.76 

.0190 

5.2 

14400 

15030 

7 

Outside  rings 

15.77 

.0200 

7.5 

16000 

15100 

8 

Outside  rings 

15.78 

.  0205 

7.5 

16800 

16030 

9 

Outside  rings 

15.76 

.0200 

7.5 

16000 

15730 

E  =  10. 5  •  106  psi 
v  =0.3 


48 


REFERENCES 


1.  W.  Nachbar  and  N.  J.  Hoff,  "On  the  Edge  Buckling  of  Axially  Compressed  Circular 
Cylindrical  Shells, "  Quart.  Appl.  Math..  Vol.  XX,  No.  3,  Oct  1962,  p.  26? 

2.  M.  Stein,  The  Influence  of  Prebuckling  Deformations  and  Stresses  on  the  Buckling 
of  Perfect  Cylinders.  NASA  TR  R-190,  1964 

3.  G.  Fischer,  "Uber  den  Einfluss  der  gelenkigen  Lagerung  auf  die  Stabilitat 
diinnwandiger  Kreiszylinderschalen  unter  Axiallast  und  Innendruck,  "  Z .  Flugwiss . . 
Jahrg.  11,  Heft  3,  Mar  1963,  p.  Ill 

4.  H.  Ohira,  "Local  Buckling  Theory  of  Axially  Compressed  Cylinders,"  Proc.  11th 
Japan  Natl.  Congress  for  Appl.  Mech.,  1961 

5.  N.  J.  Hoff,  "Low  Buckling  Stresses  of  Axially  Compressed  Circular  Cylindrical 
Shells  of  Finite  Length,"  Trans.  ASME;  J.  Appl.  Mech..  Vol.  32,  Series  E, 

No.  3,  Sep  1965,  pp.  533-541 

6.  N.  J.  Hoff  and  L.  W.  Rehfield,  "Buckling  of  Axially  Compressed  Circular 
Cylindrical  Shells  at  Stresses  Smaller  Than  the  Classical  Critical  Value,  "  Trans . 
ASME;  J.  Appl.  Mech.,  Vol.  32,  Series  E,  No.  3,  Sept  1965,  pp.  542—546 

7.  B.  O.  Almroth,  Influence  of  Edge  Conditions  on  the  Stability  of  Axially  Compressed 
Cylindrical  Shells.  NASA  CR-161,  Feb  1965 

8.  L.  H.  Donnell  and  C.  C.  Wan,  "Effect  of  Imperfections  on  Buckling  of  Thin  Cylinders 
and  Columns  Under  Axial  Compression,  "  J.  Appl.  Mech..  Vol.  17,  No.  1,  1950 

9.  W.  T.  Koiter,  "On  the  Stability  of  Elastic  Equilibrium"  (in  Dutch  with  English 
summary),  Thesis,  Delft,  H.  J.  Paris,  Amsterdam  (1945) 


49 


10.  Lockheed  Missiles  &  Space  Company,  The  Effect  of  Axisymmetric  Imperfections 
on  the  Buckling  of  Cylindrical  Shells  Under  Axial  Compression,  by  W.  T .  Koiter, 
6-90-63-86,  Sunnyvale,  Calif.,  Aug  1963 

11.  C.  D.  Babcock  and  E.  E.  Sechler,  "The  Effect  of  Initial  Imperfections  on  the 
Buckling  Stress  of  Cylindrical  Shells, "  Collected  Papers  on  Instability  of  Shell 
Structures,  NASA  TN  D-1510 


50 


NASA- Langley,  1966 


"The  aeronautical  and  space  activities  of  the  United  States  shall  be 
conducted  so  as  to  contribute  ...  to  the  expansion  of  human  knowl¬ 
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